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ABSTRACT 


One of two objectives of the present work is to 
investigate the advantage of Modal parameter extraction by 
Multidegree of Freedom (MDF) method over single degree of 
freedom (SEF) method, when the structure is highly damped 
and has closely placed modes. The MDF curve fit method used 
is a least square error minimisation technique. It has been 
tested on two structures. One was a plate held in hand and 
the other was a structure made up of four beams interconnected 
at their ends by bolts. It was found that the MDF method gave 
improved parameters compared to the parameters obtained by 
SDF method. The other objective of this work was to predict 
the response model of a coupled structure using the response 
models of the substructures. The method adopted was the 
impedance coupling method. Measured response model of the 
two substructures, a beam and a U-shaped structure has been 
used to predict the response model of the struct\are, formed 
by coupling the substructures at two co-ordinates. 
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CaA.PTER-1 

INTRODUCTION 


1 . 1 Introduction 

Measurement of the dynamic characteristics of 
materials and structures is an increasingly important part 
of the design and development of air and space crafts, 
automotive and rail road vehicles, machine tools and other 
equipment. Aided by the development of electronics and 
computers in the last two decades, dynamic modelling and 
testing has been simplified and refined considerably. 

The dynamic behavior of a mechanical structure can 
be characterised by three types of models, namely the spatial 
model, the modal model and the response model [ 1 ]. Spatial 
model characterises the physical properties, namely mass 
stiffness and damping of the system. Modal model is defined 
by the natural frequencies, damping loss factors and mode 
shape vectors. The response model is a set of responses at a 
given point due to a unit amplitude sinusoidal force applied 
to each point individually and at every frequency within a 
specified range [Ij. 

Experimental methods of obtaining the above mentioned 
models is called 'Modal testing*. The first step in Modal testing 
is the determination of the response models. This is found by 
exciting the structure, measuring the responses and determining. 
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the transfer function of the systein on the spectrum analyser. 

This constitutes the 'Hardvrare' of modal testing. The 
response model is then analysed to yield the modal and spatial 
models. This analysis constitutes the 'software' of modal 
testing. 

Analytical methods like Finite element methods, can 
also derive the mathematical models. The advantage of these 
analytical methods lies in the fact that the system character- 
istics can be determined and modified at the design stage 
itself. But the system characteristics thus predicted may 
differ from the actual performance due to complexity of 
structure, nonhomogenity of material etc. These analytical 
models can be improved by comparing it with the models obtained 
by modal testing. Modal testing can also help in monitoring 
the condition of structures and in predicting the possible 
failures of the structure before hand. 

1 . 2 Literature Review 

Literature of modal analysis dates back to 1940 's, 
when the principle of Modal testing was given by Kennedy and 
Pancu [ 2 } . The literature concerning the development of 
Modal testing is given in a recent work by Jaspal Singh [3j , 

Testing techniques can be divided into single shaker and 
multiple shaker tests. In single shaker testing, the structure is 
vibrated at one particular point and responses are measured 
at one or more points. Multiple shaker technique is used 
for large structures such as aircraft structures. A large 
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aTTiount of energy can be uniformly fed into the structures than 
the single point excitation. A number of exciters can be 
positioned judiciously at different points on the structure so 
that the modes of interest can be excited individually L 4, 5 J . 

Methods which determine the modal parameters namely, 
the natural frequency, modal damping and mode shapes, from 
the transfer function can be classified as single degree of 
freedom (SDF) methods and multidegree of freedom (MDF) 
methods. The type of method selected depends upon the 
closeness of resonances and the amount of damping. SDF methods 
extract parameters of one mode at a time. The most widely 
used method is the circle fit method [IJ • The near resonance 
points of the transfer function (Frequency response function) 
will form an arc of a circle or a complete circle, depending 
on the amount of damping. The algorithm for a least square 
error circle fit is given by Brandon and Cowley [6]. 

MDF methods may be used in both time and frequency 
domains. The MDF method given by Gaukroger, Skingle and 
Heron [7 ] is widely applicable since it can be used for 
curve fitting a large number of close resonances. The damping 
considered was of viscous type. In this method a linearised 
iterative least square procedure is followed. 

A MDF method in which many frequency response functions 
are simultaneously curve fitted, with excitation at only one 
point is given by Goyder [ 8 ] . The difference of Gaukroger 
algorithm[ 7 j with that of Goyder [8]is that the later method 
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considers each niode - ' separately, that data for each mode 
being taken from all the FRF ' s collectively. Thus a single 
mode is isolated and curve fitted for all frequency response 
functions. 

Ewins [l] has outlined a MDF method in time domain, 
called the complex exponentials method. This method is 
based on fitting a complex exponential function to the time 
impulse function. 

Most of the MDF methods require a lot of computation- A 
relatively simple method is described by Ewins and Gleeson [ 9j. 
for lightly damped structure. This method uses only the real 
part of frequency response function to determine the para- 
meters . 

The literature on coupled structures is relatively 
small. Ewins [l] has outlined a method called the impedance 
coupling method. Here the measured frequency responses of 
two separate components of a structure are used to predict 
the response of the coupled structure, formed by combining 
the two components, at one or more co-ordinates. A very 
interesting practical problem of predicting the response of 
a coupled helicopter carrier has been carried out by Ewins, 

Silva and Maleci[lO] . 

Spatial model may be easily derived from the modal 
model. A procedure to calculate the mass, stiffness matrices 
of a mechanical structure has been given by Fritzen[ll]. 
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The spatial modal derived is not alvrays unique. Interaction 
with spatial model is more meaningful if already such a 
model has been formulated and now validity or modification 
is needed. Done and Hughes [l2, 13] have analysed the 
effect of adding mass or stiffness to the existing structure 
and have established methods for obtaining the bounds within 
which the response of the modified structure will lie. 

1 . 3 Objective and Scope of Present Work. 

This work was undertaken to investigate the advantage 
of MDF over SDF method and to study the behaviour of a 
coupled structure. The MDF curve fit method was used to 
improve the modal parameters of a structure with high damping 
or closely placed modes. This iterative' m.ethod required 
initial estimates, which were obtained by SDF circle fit 
method. Two specimens, a plate held in hand and a structure 
with four beams interconnected at the ends with bolts were 
analysed by MDF curve fit method. The former was heavily 
damped while the latter had two ve3?^ closely placed modes, 
thus providing good cases to test the MDF curve fit method. 

The other purpose of this work has been to predict 
the response of a coupled structure formed by two substructures 
coupled at two coordinates. The two substructures were, a 
beam and a U-shaped structure. The method adopted was the 
impedance coupling method outlined by Ewins [l] . The predi- 
ctions were made utilizing the measured responses of the two 
substructures analysed separately. This was compared with 
the measu’~ed responses of the physically coupled structure. 




6 


The subject matter is presented in the following 
order. Chapter 2 deals with the theoretical basis of modal 
analysis. In Chapter 3, modal parameter extraction methods, 
namely SDF circle fit method and MDF curve fit method are 
described. Chapter 4 describes the impedance coupling method. 
Chapter 5 presents the instrumentation and the theory of 
signal analysis, essential for Modal testing. Results and 
discussion are given in Chapter 6. 
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Chapter - 2 

Theory of Modal Analysis 


2 . 1 Introduction 

This chapter outlines the theoretical aspects of Modal 
analysis for multidegree of freedom ( MDF) systems. The first 
section gives the basic equations of motion for a linear vibrating 
system, in time domain as well as frequency domain. The next 
section gives the response. Modal and spatial models for the 
undamped MDF system. The frequency response functions FRF ' s' 
are also defined in this section. The models for the general 
case hysteretic damping is given next. The last section gives 
the different ways of displaying the FRF data. 

2 . 2 Equations of Motion for MDF System 

The matrix formulation of the equation of motion of a 
discrete linear structure with viscous damping is 

[Mj {x(t)} +[C j{x(t)> + [Kj{x{t)} = f(t) 

( 2 . 1 ) 

where [Mj, [Cj, [K] are the mass, viscous damping and 
stiffness matrices, {x(t)} is a vector of displacements and 
I f Ctl } is the vector of the forcing functions. 
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If the forcing and displacements are assumed to be 
sinusoidal with frequency u, then {x> = {x } e^*^ and 

'k *1 

{f } = {f } e^ , and (2.1) reduces to 

[ M ] + io) [ C ] + [k] ) {x*> = {f*> , (2.2) 


The equivalent equation for harmonic motion in the 
case of hysteretic damping is 

(-to2[M] + i[K] + [Kj) {X*} = {f*> (2.3) 

The study of undamped case is discussed first because 
it is simpler than that of damped case and also aids in a better 
understanding. 

2.3 Undamped MDF System 

If in the above equations (2.2) and C2.3), [c j and [h} 
are zero then 

(.- 0)2 [ M 3 + [ K ] ) {x> = {f> (2.4) 

This is the general equation for an undamped MDF system. 
For systems under free vibration, 

^ [M ] 1 {x} 


C [ K] 


0 ) 


{ 0 } 


t2.5) 
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The solution of this eigen value problem yields 
and [{ 0 ^}^ , where is the system's natural frequency and 

{0 j-5 is its corresponding normalised mode shape vector. 

The matrix [0] possess the orthogonality property, 
which is expressed mathematically as 


[0f [M 3 [0 ] = r^i^ 


( 2 . 6 ) 


i0f [K j [0 3 = r-wM 


(2.7) 


The forced response of the undamped system is governed by 


(.[ k ]_( o ^ [ MJ 1 { x } = { f > 


C2.81 


This may be written as 


[ a Cw) 3 } 


C2.9) 


where. 


[a(wl3 = ([kJ) - 0)^ [m33~^ 


C2.10) 


Any element a , of the above matrix represents the 
response of a point j on the structure to a unit sinusoidal 
force at point Ic i.e. 


( 2.113 
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This a is known as the transfer receptance for j k and is 
jk 

called the point receptance for j = k. Also [a{(D)Jis syrnmetric 

The relation between the system's response and the force 
may be alternately expressed with respect to the velocity (v) 
and acceleration (a), as 

v. 

Mobility ¥^^( 0 ) = ^ = iu (u) (2.12) 


Accelerance a^^(a)) 


. = - 03^ a., (o) (2.13) 

fy. Jk 


The above ter^s recept.ance/ inobility and accelerance 
are generally known as frequency response functions (FRF). 

Rewriting (2.10) as, 


([k] - „^[m] ) = 


and prOTultiplying by and postmultiplying by [0] on both 

sides, one gets 

[ 01 ^ C[k] - [ M ]) [ ] =, [a Coj) 

using the orthogonality relations expressed in (2^71 and (2«,8l 

-03^)4 = [ [a((«3i [0 J 


or 


[ atol] 


2n 1 -1 r«tlT 


[0 J “ i-J i0y 


C2.14) 




of the receptance matrix [aCw)l 


is 


Thus the element a 
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a 


(oi) 


N 

2 

r=l 






(0 ) 


(2.15) 


A and 
jk 

is called the "^Modal Constant'. The equation (2.15) is then 


rewritten 


The above product 






is represented by 



(O)) 


N 

s 

r=l 



2 

(13 


(2.16) 


The Modal constants „A., , yield infonnation about 

IT J K. 

the Modal matrix [0 ] , which in turn leads to the [Mj and [K ] 
matrices by virtue of its orthogonality property given by 
(2.6) and (2.7) . 


2.4 Hysteretic Damping 

For structural or hysteretic damping, the general MDF 
equation of motion with harmonic excitation is given by (2.3). 
In absence of external excitation, the free vibration equation 
is 

C- [m] + [K] + i [H]) {X*} = {0} (2.17) 


This complex eigen value problem yields the eigenvalue 

“ic O ^ 

matrix and the eigenvector matrix [0 1 * where both the 

* 2 

matrices are complex. is the damped natural frequency 

given by 



(1 + i T^). 


(2.18) 
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where, is called the natural frequency and is the 

damping loss factor for the r mode. This [ 55 ] matrix also 
possess the orthogonality property with respect to [ M ] and 
[ K+iH ] matrices. 

For forced response the receptance matrix is now given 
by 


[ a* (oi)] = ([ k] + ifH] - [m]) ^ (2.19) 


Proceeding is the same way as in the undamped case, 

r * T 

one can express an element of [a (w) J matrix by 




N 

E 

r=l 


0 * 

jfjL 


r^k 


2 2 , . 2 
- 10 “r 


( 2 . 20 ) 


*4^ A 

Writing as the receptance is expressed as 




N 

Z 

r=l 


r^ik 


2 2 2 
_ 03 +1T^ 0)^ 


C2.21) 


2 . 5 Display of FRF 

FRF (i.e. receptance etc.) is a complex quantity and 
is frequency dependent. In order to effectively show the 
information of FRF, two types of plots are usually used 
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(a) Bode Plot: It has two parts (i) negnitude of FRF vs. 
frequency, (ii) Phase of FRF vs. frequency. 

A wide range of values of magnitude has to be shown. 

This is effectively done by using logarithmic scales for the 
modulus axis or sometimes for both the modulus axis and frequency 
axis. A peak in the magnitude plot indicates the existence of 
a mode. The frequency at which this peak occurs indicate the 
natural frequency of that mode. 

(b) Nyquist Plot: it is a graph of real part vs. imaginary 
part of FRF . It does not contain the frequency information 
explicitly. The advantage of the nyquist plot is that it 
displays the important ragion of resonance effectively. This 
plot of the resonance region will be an arc of a circle or a 
complete circle, depending on the damping. This circle called 
as nyquist circle contains the information about all the modal 
parameters namely the modal constant, damping loss factor and 
the natural frequency. This inherent information about modal 
parameters existing in the nyquist circle is effectively used 

by SDF circle fit method to estimate the modal parameters. This 
method is described in Chapter 3. 
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Chapter - 3 

Modal Parameter Extractions for a Structure 


3 . 1 Introduction 

This chapter first briefly explains the method of 
treating each mode as a single degree of freedom approxi- 
mation and then extracting the modal parameter values 
using circle fit method. The next section discusses the 
multidegree of freedom curve fit method where all the 
modes of FRF are considered together to yield the modal 
parameters. 


3 • 2 Single Degree of Freedom Analysis' with Circle Fit 
Method 

The basic assumption for the single degree of 
freedom CSDFl method is that, in the vicinity of the 
resonance, the response is dominated by that particular 
mode Cnyquist plot will be an arc of a circle) and the 
contribution of other modes may be represented by a constant 
term. Therefore the equation for receptance for hysteretic 
damping, may be written as [ i ] 

* 


a j C(j)} 


r^ik 


01 , 


7 2 

-u^+i 


+ 


* 

r jk 


C3.1) 
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The parameters ^A*^, and of the r^^ 

mode may be evaluated by following the sequence. 

Cl) Select about 10 to 15 points on either side of the 

resonance, taking care to see that the points selected 
encompass a good part of the arc of the nyquist circle. 

(2) Fit a circle to these points. The algorithm given by 
Brandon and Cowley [ 6 ] is used for this. This 
algorithm is given in Appendix A. If the quality 

of fit is better for receptance, when compared to 
that of mobility, then the type of damping included 
in the analysis is the hysteretic damping. Conver- 
sely, if the mobility gives a better fit, the type of 
damping considered is viscous damping. 

(3) A better estimate of natural frequency may be obtained 

from the point where the maximum sweep rate (de/dto^) 
occiars (Fig. 3.1). Since the points are placed at 
equal increments, finite difference is used to pinpoint 
the natural frequency with the precision of 10% of the 
frequency increment. ^ 

(4) Damping estimates are obtained by using combinations of 

2 

points below and above resonance frequency o)^, and 
the average value of all these estimates is taken as 
the damping for that mode. 

(51 The diameter, passing through of the fitted circle 
and the orientation of this diameter with the 



'Jssasi- 



Fig.3.1 Natural frequency location 
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imaginary axis, are used in the estimation of magni- 

ic 

tude of A., and its argument. The distance of the 
r j K- 

points, which is diametrically opposite to the 
natural frequency, from the origin and its orienta- 
tion with respect to the imaginary axis, is used to 

* 

calculate the term B., . 

r jk 

C6) The structures, considered in this work, are tested 
in free-free condition. Hence at low frequencies, 
the FRF exhibits mass like behaviour and at high 
frequencies, it exhibits a stiffness like behaviour. 
In order to account for these, the basic receptance 
equation may be modified by including these residuals 
as follows; 


‘jk 


+ (T- 


★ 

r^1k 


-) + 


2 M 

“ R^jk 


2 2 2 
+i 


R jk 


(3.2) 


where and ^K., are the mass and stiffness 

K J K K J K 

residuals respectively. 


The above various steps and their inherent mathe- 
matical details, has been given in a recent work [3] . 


In practice, one finds that the modes may be closely 
placed and may have very high damping. Under these condi- 
tions the effect of other mode/modes in the vicinity of a 
particular resonance is quite significant. The regeneration 
of FRF using SDF methods, may not be accurate enough. Hence 
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ths single T^od© assuniption of each resonance is not valid. 

3 .3 Multideqree Curve Fitting Method 

As explained above, the SDF methods treat data 
around each resonant frequency separately. This method 
requires simple and small amount of computation. For 
closely placed modes or for some structures with high 
damping, not enough points can be identified as belonging 
solely to a particular mode. This may necessitate other 
methods like Multi Degree Freedom (mdf) curve fitting. 

In MDF curve fitting methods, an FRF is considered 
as a single curve and the modal parameters are adjusted so 
as to achieve the best curve fit. This method, then obviously 
requires more computing power. One such method, based on 
the work by Gaukroger, Skingle and Heron [7], is briefly 
explained in the following. 


3.3.1 Details of Analysis 

, With the hysteretic damping, the receptance as a 

function of frequency is 




^ M 

« R jk 


N 

( S 
r=l 




2 2 ^. 2 
(D_' - Id +1T1., W 

r r r 


1 


+ 




(3.3) 


The parameters to be estimated in the above equation. 
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are the residuals and the modal parameters of 

th * 

the r mode - (both real and imaginary parts), o)^ 

and n^, where r = I...N 7 N = number of modes in the FRF. 

Hence the total number of parameters to be obtained are 
C4N+2) . 

* 

If Orr, is the mathematical receptance and is the 
measured receptance, then the least square error e is 
given by 

e= Z {a -a){a- a) 

erne 

over the frequency 
range of interest 

(3.4) 


_* _* * * 
where and a^, are the conjugates of and a^, respect- 
ively. This error e given by (3.4), is to be now minimised 
with respect to each of the parameters k^, where i = 1 ,...., 
(4N+2) . 

This essentially reduces to the expression given 

below (See Appendix B ) • 

[V] im + { R > = 0 (3.5) 


where the elements of matrix [P], P^j is given by 


P. . = 

3.J 




Jk^.ak^ 


.9 a 
(■ 


m 


9 a 


m 


9a 


m 


over the frequency 
range of interest 


9k 


9k. 

1 


9k, 


(3.6) 
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The elements of the vector {R}, is given by 


R - 

^i " 9k. 


over the frequency 
range of interest 


(a, - 


* 

a ) 
rn 


9 a 


9a 




m 


9k^ ^ “e 


- a 


m 


9k. 

X 


(3.7) 


Vector { 6 k} is the unknown, which is added to the 


I ‘ I I 

initial approximation k^^, k^ ^(4N+2)' 


parameters, to obtain new estimates. 

As a first step, the matrix [P] and the vector {R} 

I I I 

are formed using the initial approximations k^, k 2 ,.... ^ 4 ^ 4 - 2 * 
A good initial approximation may be taken as the modal para- 
meter values obtained from SDF, The equation (3.5) is then 

I 

solved to obtain 6 k^^ 6 k 2 , . . . values (kj^+ 6 kj^) , 

Lk 2 + 6 k 2 are then used in (3.6) and (3.7) and this 
process is repeated until an acceptably accurate solution is 


obta ined . 
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The elements of the vector {R}, is given by 



Vector { 6 k} is the unknown, which is added to the 


t ‘ I I 

initial approximation k^, k^^ ^(4N+2)' 


parameters, to obtain new estimates. 

As a first step, the matrix [p] and the vector {R} 

t I t 

are formed using the initial approximations k^, k 2 ,.-.. '^ 4 ^+ 2 * 
A good initial approximation may be taken as the modal para- 
meter values obtained from SDF, The equation (3.5) is then 

I 

solved to obtain 6 k^^ 6 k 2 , . . • values (kj^+ 6 kj^), 

Ck 2 + 6 k 2 l..... are then used in (3.6) and (3.7) and this 
process is repeated until an acceptably accurate solution is 
obta ined . 
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Chapter 4 
Coupled Structures 


4 . 1 Introduction 

This chapter explains briefly the procedure of predi- 
cting the response characteristics of a coupled structure, 
knowing the response characteristics of the substructures. 

The method used in this work is called the Impedance Coupling 
Method [ 1 ] . First this method is outlined for the coupling 

of two substructures at a single coordinate . Further on, the 
procedure for coupling at more than one co-ordinate is given. 
This section also examines, how one can include in the analysis, 
coordinates which are not involved in coupling. 

4.2 Impedance Coupling Method 

This method is also called as stiffness method [ 1 ] 

The underlying principle is given below. 

Consider two components A and B CFig. 4.1) which 
are coupled at a single co-ordinate to form the coupled 
structure C. The coupling coordinate is designated by 
numeric 1 and the coincident coupling points on substructures 
A and B are denoted by ^a and ^b, respectively. The coupling 
coordinate on the coupled structure is denoted by ^c. 

If the system A at the connection co-ordinate ^a is 

’k iost * 

excited harmonically by f e , the resulting response 


is given by 
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= a* ^ U) f* e' 
a a a 


(4.1) 


•X’ 

•1 1 


(4.2) 


Similarly for subsystem B, 


= a* ^ (CO) f* 
^b 


(4.3) 


Applying the compatability and equilibrium conditions 


at the connection co-ordinates one gets 


(4.4) 


+ fi 


(4.5) 


Substitution of (4.2) and (4.3) in (4.4) and using the result 
in (4.5), the relationship between' the receptances is 


1 1 
c c 


la la 




(4.6) 


c c 


‘1 1 
a a 


'1 1 
-^b-^fo 


C4.7) 
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where, Z denotes the impedance. 

Thus the FRF of the coupled structure C can be obtained 
in terms of the FRF properties of the substructures analysed 
separately. 


The above analysis can be extended to the systems conne- 
cted at more than one co-ordinate. The co-ordinates which are 
not involved in coupling may also be included. The whole process 
involves the partition of the subsystem impedance matrices 
into submatrices which are involved in coupling and those which 
are not and then adding together only those submatrices which 
are involved in coupling. 

Therefore, if 3 [z^{w)l are the recep- 

tance and impedance matrices, respectively, of subsystem A 
and [ag(tol] and [ZgC‘'^)]are the corresponding matrices of 
subsystem B, and if the coupling coordinates are collectively 
designed by q, the coordinates not involved in coupling by 
p for subsystem A, then the receptance matrix is partitioned 
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[z (U)] = [a,(“)r^' 


PaPa 


^aPa 


P <3 
^a^a 




(4.9) 


Similarly the impedance matrix for the subsystem B, 
is given by 


[ Z* Co))] 


* -1 
[ 0!g(«3) ] 


^b^b 


^b^b 


^b^b 




C4.10) 


where, r denotes the coordinates on subsystem B which are not 
involved in coupling. 


As before, applying the compatibility and equilibrium 

* 

3 

c 

structure can be obtained in the form. 


condition, the impedance matrix [z (.w)3 for the coupled 


[Z*(a3)] = [Z*(a))] + [Z*(a))] 


P P 


‘^aPa 


Vb 


Vb 


P q 
a^a 




★ * 

z +z 

^a^a ^b^b 


(4. HI 
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The receptance niatrix combined 

★ 

systern is obtained by the inverse of [ (o))] as 

[a* (o)) ] = [ (w) ] (4.12) 

From modal testing point of view, the impedances 
required for the substructures may be evaluated by any one 
of the options mentioned below. 

Model . 

(i) Using the modal^and evaluating the receptance matrix 

using the formula 

* *_ 2.i-lr_»*TT 

[ a (to) J = [{ZT J )'4 [0 J 

(4.13) 

where 

X* = (0^ (1 + iTij.) 

(ii) Calculating the inverse of the receptance FRF matrix 
obtained by direct Ttieasureinent • 

(iii) Calculating the inverse of the regenerated receptance 
FRF matrix. 

The results presented in section 6 are by using the 
last method of utilizing the regenerated FRFs. 
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Chapter - 5 

InstruniePtation and Signal Analysis 


5 . 1 Introduction 

Experimental modal analysis involves exciting a 
system, measuring the excitation and response time signals 
and then operating on these time signals to determine the 
system transfer function. The instrumentation used in this 
work is shown by the block diagram in Hg. 5.1. The impact 
hammer, the accelerometer and the preconditioning amplifier 
are first described briefly in the followings. The signal 
analysis is later described in section 5.2. 

Excitation by Impact hammer is one of the commonly 
used excitation technique . The impact excitation does not 
not require elaborate fixture and it is ideal for use in 
restricted spaces where an exciter would not fit. Structures 
of masses in the range of approximately 2 kg to 1000 kg can 
be excited by impact -hammer . The equipment consists of an 
impactor with different types of tips and heads, which serve 
to extend frequency and force level ranges of excitation. 

Just below the tip of the impactor is a force transducer, 
which detects the magnitude of the force felt by the impactor. 
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Accelerometers are piezoelectric transducers 
used to measure the response (acceleration) of a system. 

The piezoelectric crystal attached to a seismic mass, senses 
the inertial force of the seismic mass. This force gener- 
ates an electric charge across the end faces of the crystal, 
which is proportional to the acceleration felt by the acce- 
lerometer. Accelerometers are preferred over other types 
of pickups, namely the velocity and displacement pickups, 
because accelerometers are much smaller in size and can 
cover a wider frequency range. 

Amplifiers are used to boost the very small 
electrical charge that is generated by the piezoelectric 
transducers into a signal strong enough to be measured 
by the analyser. There are two types of amplifiers namely 
charge amplifier and voltage amplifiers. Charge amplifiers 
are preferred over voltage amplifier, mainly due to the 
reason that the sensitivity and gain are not affected by 
the length and properties of the connecting cable in case 
of charge amplifiers, but is adversely affected in case 
'of voltage amplifiers. 



Charge 

Amplifier 

2626 



Impact Hammer 
8202 



Test piece 


Accelerometer 

4369 


Charge Amplifier 
2635 

Signal Analyser 
SM-2100B 


□ 


Channel i 
Channel 2 


P.C.-XT 


IEEE-488 
(G.P.I B.) 


HIOKI Recorder OR 0 ^ » 5 oo 

8801 X - Y J^ecol^4e^ 


Block Diagram of Experimental Setup 
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The spectrum analyser receives the two time 
signals namely, the excitation and response signals. The 
analyser operates on these time signals, to determine the 
system transfer functions or frequency response function. 
The specifications of the instruments used are indicated 
in Fig. 5.1. 


5.2 Signal Analysis 

5.2.1 Fourier Transforms 

The first step in modal testing is to relate the 
input (excitation f (t) ) and the output (response x(t)) of 
the system in the frequency domain. This is done by performing 
dual channel FFT analysis on a spectrum analyser. 

If f (t) and x(t) are the input and output of a 
system in time domain, then the function relating these two 
in time domain is called the Impulse Response Function h(t). 


This relationship is given by 
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x(t) 


-f oe 


f 


h (t ) f (t, -T ) dr 


(5.1) 


In frequency domain this relation is 


X*{co) = H*(a)) . F*(a)) (5.2) 

The steps involved in arriving at (5.2) from (5.1) are given in 
Appendix c 

ic * 

In (5.2), X (u) and F (a>) are the fourier transforms of 

★ 

xCt) and f (t) respectively and H (lo) is the frequency response 
function of the system. 


then 


If x(t) and f (t) satisfy the Dirichlets' condition. 


and 


* 

X (u) 


1 _ 

2Tr 


+ CO 


/ x(t) 



F* CojI = 


2 Tt • 


; f(t) e"^“^ dt 


(5.3) 


(5.4) 


5.2.2 Autospectra, Cross Spectra, Transfer Function and 
Coherence 

5. 2. 2.1 Auto Spectra 

Auto spectra of x(t) and f (t) are defined as 


S (.0)1 

XX ^ 


ic ic 

X Co)l - X (tol 
% 


^ff ~ 


F* CojI . F* (u) 


(.5.5) 
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where x im) and F (w) are the complex conjugates of X (03) 

'\j ^ 

•k 

and F (u) respectively. The autospectrum is a real quantity 
and it indicates the distribution of power in a signal as a 
function of frequency. 

5. 2. 2. 2 Cross-Spectrum 

It is the basic function which correlates the input 
and output signals. The cross -spectrum between the two signals 
x(tl and f (t) is given by 



li 

e 

• 

X* (to) 


k 

k 

* 

(5.6) 


= F (to) . 

X (to) 

* 

* 



S 

xf 

= 

-vfx 



5 . 2 . 2. 3 Transfer Function 



This 

function relates the input 

and output of the 


system. Thus from (5.2), the transfer function as a function 
of frequency is 

(5.7) 

F (o)) 

Alternately, the transfer function is defined in terms of 
the autospectrum and the cross spectrum as, 

H Oo) = ^ (w) 


(5.8). 
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or 

H*(o 3) = / S*j{a)) (5.9) 


The division in (5.8) is by a real quantity which 
is rnore convenient on the spectrum analyser. 


5. 2. 2. 4 Coherence 


For exact measurements, (5.8) and (5.9) should yield 

2 

identical results. Based on this a coherence function y (oj) 
is determined as 




• * f 2 


(5.10) 


The coherence is ideally equal to one for linear, 
deterministic, noise free system and is between zero and one 
otherwise. Thus, the coherence forms a good measure of faith- 
fulness of the experimental data. Later in Chapter 6, the 
coherence is utilized to justify the measurement of acceler- 
ance rather than receptance. 

5.2.3 Discrete Fotirier Transform DFT 

With the advent of digital computers, the focus is 
entirely on discretization of analog signals. In evaluating 
DFT, both the time and frequency domain data are in discrete 


form 
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The tirne signal x(t) recorded for a sample length T, 
is discretised to obtain N discrete evenly spaced values. 

The sampling rate, u is given by 

o 


(s)g = 2 itN/T 


(5.11) 


* 

DFT, X (to) , of this time signal, may also have N 

discrete values. If so, then the frequency range is zero 

to w , with a frequency resolution of Aai, where 
s 




tOg/N. The forward transform will have the form 


* 

X (k) 


i 

N 


N-1 

Z x(n) 
n=0 


-i 


e 


2TTkn 

N 


k = 0, 1, . . .N-1 

(5.12) 


and the inverse transform takes the form 


% ^nl 


N-1 * 

z X (k) 
k=0 


. 2TTkn 

XT 


n = 0, 1. ..N-1 

(5. .13) 


Equation (5.12) is expressed as, 

{X*> = I [A*3 {x> (5.14) 

'fc 

where {X } is a vector representing the N complex frequency 
components, [ A } is a square matrix of unit vectors (order NxN) 
and {x I is a vector representing the N, discrete time signal 


values 
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Equation (5.12) indicates that obtaining N frequency 
components from N time samples requires N complex multi- 
plications. Using fast Fourier transform (FFT) algorithm given 
by Cooley and Tukey [l4/ same result can be obtained with 

a greatly reduced number of complex multiplications (of the 
order of N Log^) . 

5.2.4 Aliasing 

This problem is associated with digital spectral 

analysis. If the sampling rate, co , is too slow, then the 

s 

high frequencies in the original time signal may be misinter- 
preted as low frequencies. This form of error is known as 
aliasing. This problem of aliasing is illustrated in Fig. 5.2. 
Figure 5.2 (al shows the case in which the sampling rate u 

s 

is higher than the frequency of the sinusoid. Hence on 
performing DFT, the actual frequency is correctly recognised. 
For the case shown in Fig. 5.2 (b), the sampling rate, oo-,/ 
is lower than the frequency of the sinusoid. As can be seen 
from the same figure, the sampling points are now situated in 
such a way that on performing DFT, a lower frequency is 
recognised instead of the actual frequency. 

The relationship between the sampling rate and 
the frequency content of the time signal is expressed in 
Shannon's sampling theorem [I 4 ] ^which states that a sampled 
Cdiscretisedl time signal must not contain components of 
frequencies above half the sampling rate C called the nyquist 
frequencyl. The maximum frequency contents of the time signal, 
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“max 


is. 


“max 


»^/2 . 


(5.15) 


and the resolution of the frequency spectrum, Aj) is 


= (i) „/N 

s 


(5.16) 


It should be noted here, that with this restriction 
of Tnax' number of discrete frequency points is N/2. 

Equations (5.12), (5.13), then get modified as 


^ N_i 

X (k) = i I x(n) e ^ , k = 0, 1.. . cf - 1) 

n=0 ^ 


* 

X Cnl 


N/2 - 1 

2 xCt) 

k=0 


i 


e 


2Trkn 

N 


X = 0, 1. . . (N-1) 


(5.17) 


In the frequency domain, the error caused due to 
alising is depicted in Fig. 5.3(a). If the time signal does 
have contents of high frequency then the signal is passed through 
an anti -alising filter to cut off the higher frequency 
signals. This cutoff frequency is decided by the frequency 
of sampling rate, One such filter and the filtered signal 

spectrum is shown in Fig. 5.3 Cb). 

5.2.5 Leakage 

The time signal' is recorded for a finite length of 
time history and the signal is assumed to be periodic in 
this time length. .'f f Is signal is not fully periodic, 
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Alias Distortion of Spectrum by DFT 

F1G.5-3(q) 





Anti-Aliasing Filter Proccess 

FIG. 5-3(b) 
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then the problem of leakage appears. This can best be 
illustrated by the following example. Two samples of the 
same signal are recorded at slightly different record ' 
lengths. The signal shown in Fig. 5.4(a) is fully periodic 
^nd hence the DFT indicates the single frequency shown in 
the same figure. The signal shown in Fig. 5.4(b) is not fv 
fully periodic and the DFT does not show a single frequency 
as in Fig. 5.4(a). This indicates that the energy has leaked 
into a number of the spectral lines close to the true frequency. 
When the signal frequencies are low, this leakage problem 
can be greatly disturbing. The problem of leakage can be 
solved by proper 'windowing', described in the next subsection. 

5.2.6 Windowing 

Windowing involves the weighting or imposition of a 
prescribed profile on the time signal x(t) prior to performing 
the Fourier transform. The weighting is denoted by W(t) . The 

I 

analysed signal x (t) is given as 

X' (t) = x(t) . W(t) (5.17) 

Some common types of windows used are shown in 
Fig. 5.5. The Iwatsu spectrum analyser offers rectangular 
and Hanning windows and there are described briefly in the 
following. 



Wn or n 

Leakage of Spectrum 
FIG. 5-4 (b) 





Fig. 5.5 ! 


Some Common Types of Windows 
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(i) Rectangular Window 

The rectangular windowing or weighting, also called 
tlat or Boxcar weighting, in actually no weighting, all 
on the finite tiine record (T) . it is defined as 


W(t) = 1 for 0 < t < T 

W(t)’ = 0 elsewhere. 


(5.18) 


The practical use of the rectangular window is for 
analysing transients wdth shorter durations then the record 
length T. Rectangular windows are easily the first choice for 
all types of signals. 


(ii) Hanning Window 

Hanning window is a sum of rectangular window and period 
of a cosine of equal amplitude. Hanning weighting is defined 
as 

W(t) = 1 - cos 2n:t/T for 0 < t <_ T 

C5.19) 

W(t) = 0 otherwise. 

Hanning windows are typically used for continuous signals 
such as produced by steady periodic or random vibration. It 

0 tf©ctively the problems of leahage. This is illustrated 
by using the banning windowing for the example given in Fig, 5.4. 
The improvement gained by using banning windowing is evident in 
Fig. 5.6. 
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FIG. 5-6 
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5.2.7 Averaging Modes 

Averaging of a number of samples of time signal 
are done usually for processing random signals. When 
analysing random vibration, it is the estimates of the 
spectral densities and correlation functions, which 
actually characterise the type of signal. Three modes of 
averaging are usually provided in the spectrum analyser. 

(a) Sum Modes This operation is executed by assuming 
the weight of each signal as 1/N. The noise of the 
averaging signal is decreased with the increase of 
the number of samples averaged. 

(b) Exponential Mode: This operation is executed by an 
exponential weighting of the samples of time signals. 

In other words, the starting samples are given more 
weighting than the latter samples. 

Cc) Peak Mode: This operation compares the size of each 

sample of time signal to select the sample with 
maximum value. 

5 .2 »8 Iwatsu Spectrum Analyser and the Options Selected 
The front panel of Iwatsu Spectrum Analyser offers 
a large class of analysing functions by operating function 
selection keys. function used in this work was coherence 

function. The use of this function automatically also 
calculates autopower, crosspower ancl the transfer functions, 
which are registered in separate blocks by the inbuilt 
functions. These autopower functions could also have been 
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5.7 (a) 


Time history of impulse excitation signal 
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Fig. 5.7b 


Time history of response signal ( 2.J<:) 
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Fig. 5.7b ; Time history of response signal (4K) 
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Fig. 5.8 : 


FRF obtained from response signals of 
different record lengths 
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calculated separately by individual processing keys. 

The time signal is recorded by selecting the data 

length (i.e., N) and the upper limit of analysing frequency. 

The latter choice automatically fixes the pass band frequency 

of the antialiasing filter. The sampling rate, w , is 

s 

decided by (5.16). 

As explained in Section 5.2.6, a proper choice of 
window function is needed. This is to be guided by the type 
of signal to be recorded. The impact hammer produces an 
impulse signal of very short duration, the signal being 
produced at the beginning Instants, as shown in Fig. 5.7(a). 

This calls for an rectangular window as discussed in Section 

5.2.6. The two channels of the analyzer are guided by the 

same window. A typical response signal with data length 

as 2K and upper frequency range as 2 KHz is shown in Fig. 5. 7(b) 

and it is clear that wdth the record time of 250 msec, the 

signal has not decayed completely. In order to have faith in 

in the results, the response signal was also recorded over 

500 msec with data length as 4K, and upper free range as 

2 KKz. The resulting vrf's using these two signals 

are compared in Fig. 5.8 and it is evident that the record 

time of 250 msec is sufficient to yield the correct information. 

The extraction of modal parameters from the recorded 
data needs post processing of the FRF's. For this post 
processing the FRF data registered in spectrum analyser is 
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transferred to the IBM PC/XT via GPIB (IEEE - 488) bus. 

The results are later transferred back to signal analyzer 
for the display and recording purposes. The data transfer 
conrimands are executed from the PC/XT. 
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C-hapter fi 
Results and Discussion 


6. 1 Introduction 

A previous work [i ] reported results using Modal analysis 
through a SDF method, namely the circle fit method. This chapter 
first discusses results for a hanging plate held with hand to 
illustrate the inadequacy of the circle fit method when damping is 
high, eventhough the modes are well separated. An MDF curve fit 
method is used to modify the modal parameter values and the MDF 
method is found to generate more accurate results. The next section 
discusses the problem of a structure with last two of the five modei 
being closely spaced and the MDF method is again shown to yield 
better results for the modal parameters. 

The last section presents the results for a coupled struct- 
ure. Here the two substructxires are separately analyzed for modal 
parameters and the FRF‘s for receptance of the coupled structure 
is generated using the values of modal parameters of the substruct- 
xares. This is checked against the experimentally obtained FRF for 
the same coupled structure. 

The FRF's displayed in the chapter are for accelerances 
but the modal parameter extraction are done from the receptances. 
The accelerance is measured instead of the required receptances* 
Receptances are calculated by (2.13). This additional computation 
required is due to the fact that the double integration process 
in charge amplifier results in the usage of high gain in amplifier 
which introduce noise . This can be clearly seen by observing 
coherence for both acceleration and receptance meastarement [sj. 
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A case where the SDF fit method gives inaccurate 
estimates of the Modal parameters has been studied. The 
structure was a plate which was held lightly by hand to 
induce high damping. Upto a frequency range of 500 Hz, 
this structure exhibited two modes, of which the first 
is highly damped. A regenerated FRF obtained by SDF circle 
fit method superimposed on the measured FRF is shown in 
Fig. 6.1, where the experimental curve is shown by the 
jagged natured plot. From this figure it is evident that 
the SDF circle fit m- 'd has given a bad fit for the first 
mode. With MDF curve fit method, described in Chapter 3, 
the regenerated FRF has been improved. This structure was 
used mainly to illustrate the improvement of modal para- 
meters, mainly the damping ratios. So the residual terms are 
neglected and thus eight parameters need to be considered 
in C3.51. The initial estimates required for MDF curve fit 
method has been taken from the parameter values provided by 
SDF circle fit method. 


The number of iterations and error after each iter- 
ation is given in Table 6.1. For this case, the results 
converge after only three iterations. The Initial estimate 
of parameters taken from SDF circle fit method along with 
their refined values from MDF curve fit method after 
5 iterations are shown in Table 6.2 The regenerated FRF 
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Fig. 6.1 : Measured and re erated (SDF ^method) FRF for 
plate held in hand. 



Fig. 6.2 s Measured and regenerated (MOF method) FR- 
foiT Tia.iia. 
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Table 6.1 : Error after each iteration of MDF curve 

fit method for a plate held in hand 


Error after SDF circle fit method = 72.942 


Iteration No, 

Error after each 
iteration 

1 

11.094 

2 

7.458 

3 

7.413 

4 

7.408 

5 

7.407 


Time taken for each iteration ^ 100 seconds 
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Table 6.2 : Parameters obtained by SDF circle fit method 

compared with improved values obtained by 
MDF curve fit method for the plate held in 
hand 


Parameters 


SDF circle fit 
method 


MDF curve fit 
method 


n 

1 

“1 

"k 

^l^jk^real 

k 

imag 

n2 

“2 

k 

^2^jk^real 

k 

^2^jk^imag 


0.042665 

217 .8 
0.166144 
n. 008244 
0.009244 
439.3 
-0.091384 
-0.011683 


0.0721406 

278.11051 

0.279986 

0.052724 

0.0C92857 

439.35958 

-0.09195 

-0.008507 
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obtained from these improved values superimposed on the 

measured FRF is shown in Fig. 6.2. Comparing Fig. 6.2 

with Fig. 6.1, it is evident that the MDF curve fit method 

yields significantly better results. The improvement of 

the first modal damping constant as well as of the 

* 

modal constant, noticable from Table 6.2. 

6 . 3 Numerical Procedures and Results for a Jointed 
Structure 

6,3.1 Jointed Structure Configuration 

The structure is made up of foxir beams interconnected 

at the ends as shown in Fjg. 6.3. This structure offers an 

advantage that by changing the included angles betw'een the 

limbs, the position of the modes along the frequency axis can 

o 

be shifted. This structure with included angles as 68 and 
112°, shown in Fig. 6.3, was analyzed in an earlier work [3 ]. 

When the included angles are not near 90°, the modes 
are well separated for a frequency range upto 500 Hz. This 
can be clearly seen from the FRF of Fig. 6.4, which has been 
reproduced here from [ 3] , and exhibits five well separated 
modes. The modeshapes of these five modes are also 
reproduced in Fig. 6.5. 

In the present work, the same structure has been 
studied when the included between the limbs are very close to 
90°. The structure with this configuration is shown in 
Fig. 6 . 6 . A typi-cal FRF a* 4 , obtained by measuring the 
response at 1 with impact excitation at 4, is shown in 
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* 

: FRF a^jL 


Fig. 6.4 


of Rhombus shaped Jointed Structure 
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Mode 3 



Mode 4 



Mode 5 


Fig. C-5(h) 


MODH S’-^rvPES OF R-’OjV3US SH.-.PED STRUCTURE 
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Fig. 6.7- The fourth and fifth modes are now seen to 

be quite close to each other. The frequency resolution is 
1 Hz and even with this resolution it is difficult to take 
a sizeable number of data for SDF circle fit method, without 
zooming facility in the signal analyzer. The modal para- 
meters were calculated using SDF circle fit method and 
were used to regenerate the FRF. This regenerated FRF, 
superimposed on the experimentally obtained FRF is shown in 
Fig. 6.8. It is clear from this figure that the SDF 

circle fit method does not result in a good regeneration in 
the region surrounding the fourth and fifth modes, whereas 
upto the third mode the regeneration is acceptable. The 
sinale mode assumption at these last two modes is not really 
good since each of these two modes will have a significant 
effect on each other. Moreover, the number of points required 
for a circle fit, available in this case is very small, which 
leads to an erroneous circle fit. This is turn leads to 
erroneous estimation of modal parameters. This configuration 
indeed gives us a good case for testing the MDF curve fit 
method . 

6,3.2 Improvement of Modal P arameters by Curve Fi ^ 

Method 

The MDF curve fit method, described in Chapter 3, 
was applied to the jointed structure with configuration 
described in 6.3.1, the structure having exhibited closely 
placed fourth and fifth modes. 
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Fig. 6.7 
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Equation (3.11) Is reproduced here ae 

[P 3 {6k} = (R} 

where the individual elements of [p J matrix, P. . and 
vector {R}, are given by (3.6) and (3.7) respectively. 

These values are calculated over the frequency range using 
the previous estimates of the modal parameters. From (6.1), 
{ 6k } is then obtained; which represents the incremental 
change in the modal parameter value vector {k} . These 
increments { 6k > , are added to the previous estimates to 
obtain the new estimates, i.e. 


I 


f 

Previous j 

1 

New estimate 

J estimate of 

+ { Sk) = < 

of 

parameters j 

Mr 


^ parameters 


} 

( 6 . 2 ) 


It was observed that on regeneration of FRF using these 
new estimates, the error given by (3.4), increased with each 
iteration drastically indicating the divergence of the 
solution. 

It was also observed that on scaling these increments 
{ 6k} by a factor e , where e •<; 1, the solutions converged 
monotonically. Therefore the new vector {k} is obtained by 




Previous ) 

, New estimate 

estimate of V t e {6k} = , 

of parameters 

parameters j ! 

V 


(.6.31 
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The value of e 'tiay not be the savie for each 
iteration. It is chosen by trial and error, heepincr in view 
that the error obtained by using the new parameters is lesser 
than the error obtained by using the previous estimates. 

The explanation for the divergence for e = 1 and 
convergence for e < 1 is due to the following reason. The 
solution for these FRF ' s with closely placed modes has a 
very small stable solution region around the actual solution. 
If the solution crosses this stable region and reaches the 
unstable region, the method diverges. 

6.3.2. 1 MDF Curve Fit Method Applied for Full Frequency 
Range and all the Parameters 

The MDF curve fit method with scaling factor £ was 

•k 

used to obtain the modal parameters for the FRF shown 

in Fig. 6.7. The initial estimate vector v/as taken as 

the parameter values obtained from SDF circle fit method. 

Table 6.3 lists the scaling factor e and the error e after 
each iteration. Table 6.4 compares the modal parameters 
as obtained from SDF circle fit method and as obtained from 
the MDF curve fit method after 17 iterations. It is evident 
that the parameters associated with only the fourth and fifth 
modes get refined. The regenerated FRF using the MDF curve 
fit results is shown in Fig. 6,9. On comparing this with 
the regenerated FRF of Fig. 6.8, the improvement near the 
last two modes is evident. This is better illustrated if 
the nyquist plots obtained from SDF and MDF results are 
compared in Fig. 6.10. These nyquist plots are drawn only 
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Table 6.3 : Error and scaiir, ^ ^ 

of MDF curve factor for each iteration 

FRF a* iTiethod for jointed structure 

14 

Error after sDr circle fit method = 323.495 


Iteration 

number 

Scaling 
factor e 

Error after the 
iteration 

1 

O.Ol 

307.933 

2 

0.02 

282.319 

3 

0.03 

253.858 

4 

0.04 

226.611 

5 

0.05 

201.634 

6 

0.06 

178.612 

7 

0.07 

156.992 

8 

0.08 

136.523 

9 

o.io 

115.327 

10 

0.10 

98.139 ' 

11 

0.10 

84.139 

12 

0.20 

62.581 

13 

0.30 

42.618 

14 

0.30 

32.299 

15 

0.30 

26.922 

16 

0.30 

24.088 

17 

0*30 

22.782 

Time taken for 

each xteratW^: on 
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Table 6.4 : 

SDF versus MDF for 

Jointed Structure FRF at , 

14 

Parameters 

From SDF circle 

From MDF curve fit 


fit method 

method after 17 iter- 
ation 

Mass term 

0.26173388 

0.28799294 

n 

1 

0.01412377 

0.01474284 

CD 

1 

103.4 

103.374934 

ic 

4^14 ^real 

^0.08636302 

-0.09076631 

★ 

^1^14^imag 

0.02284896 

-0.00158237 

^2 

0.00067156 

0.00067967 

ui2 

207.5 

207.55378 

•k 

^2^14^ real 

-0.2558857 

-0.255804 

k 

^2*‘14^imag 

-0.050539 

-0.050442 

^3 

0.00163598 

0.00231796 

“3 

288.7 

288.7808 

k 

^3^24^ real 

0.03305854 

0.03002156 

k 

^3^14^ imag 

0.01132084 

0.01052005 

^4 

0.00906141 

0.00956932 

“4 

365.8 

365.93098 

★ 

^4'^l4^real 

-0.143115 

-0.1383204 

k 

^^I4^imag 

-0.024972 

-0.043329 


0.0087456 

0.0095888 

“5 

374.7 

374.96568 

k 

^^14^reak 

0.071151 

0.0500657 

k 

^5^14^ imag 

0.045904 

0.0137109 


70 


from 300 Hz to 450 Hz so as to signify the results around 
fourth and fifth modes. 

taken for each iteration was about 20 minutes 
on the PC/AT machine with 8087 chip. This is an enormous 
computation time considering the fact that MDF method may need 
many iterations to converge. The next subsections discuss two 
alternatives to cut down this computation time. 

6. 3.2.2 Reducing the Number of Parameters to be Improved 

As mentioned earlier in 6.3.2. 1 and illustrated 
in Fig. 6.8, the SDF circle fit method gives a satisfactory 
estimates for the first three modes. The major error is 
indeed from the estimation of modal parameters of the fourth 
and fifth modes. So the MDF curve fit method has now been 
applied for the same FRF 3^4' improve only the p)arameters 
03 r' real part of and imaginary parts of where 

r = 4 and 5. The frequency range covered is the full range i.e. 
upto 500 Hz. Since the number of parameters to be improved are now 

only 8, the time required is brought down to only about 3^ 
minutes per iteration. The scaling factor e and the error 
after each iteration is given in Table 6.5. The regenerated 
FRF superimposed on the measiared FRF is given in Fig. 6.11 
and the regeneration is observed to be quite faithful. The 
nyquist plots also were seem to match, ^falile compares 

the parameters obtained by MDF curve fit method with the 
parameters of SDF circle fit method. 
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Table 6.5 


scaling factor e for each iteration 
of MDF curve fit method of section 6. 3.2.2 

for J ointed structure FRF 


Iteration 

No. 

Scaling factor 

6 

Error after 
each iteration 

1 

0.01 

310.189 

2 

0.03 

279.37 

3 

0.04 

254.304 

4 

0.08 

224.952 

5 

0.08 

207.889 

6 

0.12 

188.944 

7 

0.08 

179.203 

8 

0,08 

170.921 

9 

0.08 

163.863 

IQ 

0,1 

156.425 

11 

0.1 

150.326 

12 

0.1 

145.316 

13 

0,1 

141.195 

14 

0.12 

137.167 

15 

0.12 

133.976 

16 

0.12 

131.44 

17 

0.20 

128.23 

18 

0.30 

125.147 

19 

0.3 

123.46,4 

20 

0.1 

123.122 


Time taken for each iteration minutes. 
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Table 6.6 : SDF versus MDF of section 6. 3. 2. 2 for 

•k 

FRF of the Jointed structure 


Parameters 

SDF circle fit 

MDF curve fit 


method 

method 


0.00906141 0.0094403 

0)4 365.8 365.942627 

^4^l4^real -0.143115 -0.13690983 

(Atj. -0.024972 -0.04398023 

Y 14 inag 

n_ 0.0087456 0.00946401 

0 

0 )^ 374.7 374.99396136- 

b 

0.071151 0.04734567 

* 

^5^14^ imag 


0.045904 


0.01151493 
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Fia. 6.11 : Measured and regenerated (MDF jiethod - restricting 

number of parameters) of FRF a^^ of Jointed structure 



Measured and regenerated {MDF method - restricting 
both frequency range and the number of parameters! 
of FRF a^^ of Jointed structure 


Fig. 6.12 : 



6*3.2. 3 Restri cting the Frequency Range and the Nuniber of 
Paratneters to be Improved 

The contribution froixi the parameters of the fourth 
and the fifth modes to the frequencies away from these 
modes are found to be extremely small. Hence the frequency 
range covered has been restricted to that around the fourth 

and fifth modes, i.e. from 325 Hz to 410 Hz. 

★ 

The results for the same FRF shown in 

Fig. 6.12 and Tables 6.7 and 6.8. As can be seen from 
Table 6.7, the error has been reduced to almost the same 
level as in the Table 6.5, though the frequency range now has 
been reduced. 

The advantage of this restriction of frequency range 
is evident from the fact that the time taTcen for each vibration 
is now only 40 seconds compared to the iteration time of about 
3 J 5 minutes for the results of Table 6.5. 

6.3.3 Mode Shapes for the Jointed Structure 

For obtaining the mode shapes for the jointed 

structure, with included angles very close to 90^, a total 

of 13 FRF were measured, one point FRF and 12 transfer FRF‘s, 

* * * . * 

Out of the 13 FRFs, four FRF‘s a^^, a^^, and a^g respon- 

ded well for SDF circle fit method which is evident from the 
Fig. 6.13. The fourth and fifth modes appear to contribute 
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^^<3 scaling e for each iteration 
? ^'^fve fit method of section 6. 3. 2. 3 

tor RF Jointed structure 

Error after SDF circle fit method: 323.495 


Iteration 

number 


^"^aiing Error after 

factor e the iteration 


1 

O.Ol 

310.186 

2 

0.03 

279, 368 

3 

0.04 

254.302 

4 

0.08 

224.951 

5 

0.08 

207.888 

6 

0.12 

188.943 

7 

0.08 

179.202 

8 

0.08 

170.92 

9 

0.08 

163.863 

10 

0.10 

156.425 

11 

O 

0 

O 

150.326 

12 

0.12 

144.384 

13 

0.12 

139.691 

14 

0.12 

135.977 

15 

o 

0 

a 

00 

133.969 

16 

0.08 

132.245 

17 

0.06 

131.121 

18 

0.06 

130.116 

19 

0,06 

129,216 

20 

0.12 

127.656 

21 

0.20 

125.669 

22 

0.3 

123.746 

23 

0.4 

122.406 

24 

0.6 

121.695 

25 

0.6 

121.389 


Time taken for each iteration 40 seconds 
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Table 6.8 : SDF versus MDF of section 6. 3. 2. 3 

for FRF of the jointed structure 


Parameters 


After SDF circle After MDF curve 

fit method fit method 


^4 

“4 

ic 

^4'^14^ real 
* 

^4^14^ imag 

“5 

•k 

^^14^real 

* 

^^14^ imag 


0.00906141 

365.8 

- 0.143115 

- 0.024972 

0.0087456 

374.7 

0.071151 

0.045904 


0.00943072 

365.93578896 

- 0.13652802 

- 0.04460380 

0.00936464 

374.97755729 

0.04671169 

0.01176846 
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(a) FRF a*j 



' (b) FRF a*7 

rig. 6.13 = Measured and regenerated (SDF a,ethodl FRFs of Jointed 
Structure 
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(c) FRF a*g 



(d) FRF a* 3 



Figure 6.-13 
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equally for these FRFs, thus providing sufficient number of 
points to do the circle fit. The other FRFs were all impro- 
ved over the region covering the fourth and fifth modes by 
the MDF curve fit method described earlier in section 
6. 2.2.3. The SDF circle fit method did not yield good esti- 
mates for the modal parameters of fourth and fifth modes. 
Thus the modal parameters corresponding to only the fourth 
and fifth modes were refined using the MDF curve fit method 
and the resulting regenerated FRF*s are shown in Fig. 6.14. 


The modal constants A*, do not explicitly yield the 

r JK 

* 

elements of the normalised eigen matrix [0 J . In order to 
extract the individual elements of this matrix , one needs 

at least one point FRF and other transfer FRFs. These transfer 
FRF's should have a common excitation point or response point 

•k 

as that for the point FRF. The modal constant is given 

by 


•k k k 

r^jk r^j • r^^k 


(6.4) 



From (6.5) above, the element can be calculated. 

This corresponds to the k^^ row of [0 j matrix. For other 
transfer FRF, say a*. , with same excitation (or response) 

point, the row of the [0^^ matrix can be evaluated by 
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P 


(a) FRF a * 2 




Figure 6.14 : Measured and regenerated (MDF method) FRF's 

of Jointed structure 
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* 

1-12 




FRF 


> a 


* 

1-13 


Figure 6.14 
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using the relation 

r^j ~ r^jk ^ r^k (6.6) 

using (6.6) above, all the elenients of [0^can be calculated. 

These elements of the eigen matrix are calculated from 
the previously obtained modal parameters using SDF circle fit 
method combined with MDF cvurve fit method. This eigen matrix 
is given in Tat)le 6.9. The average of natural frequencies 
and damping obtained for the thirteen FRFs for the five modes 
are given in Table 6.10. 

The eigen matrix is complex indicating that the 
motion is not ‘Synchronous* in a damped system. Since the 
imaginary parts are relatively small compared to the real 
part, the mode shapes are drawn using only the real parts of 
[0 ] ponsidering a synchronous motion. These mode shapes are 
shown in Fig. 6.15. 

On comparing these mode shapes with that of the mode 
shapes obtained when the included angles are 68° and 112° 

(Pig. 6.5), one finds that all the five mode shapes are 
similar for both cases. 

6.4 Coupling of Structures by Impedance Coupling Method 

The method outlined in Chapter 4 has been applied to 
the coupling of two substructures - a beam and a 'U* shaped 
structure. These two substructures and the resulting coupled 
structure are shown in Fig. 6.16. The coupling of the two 



Table 6. 
Mode 1 


Mode 2 


5 : Mode Shape Vectors for 


Real Part 
0.3135 
-0.0979 
-0.4364 
-0.2758 
-0.1816 
0.1193 
0.3502 
0.1240 
-0.0233 
-0.2767 
0.0616 
0.2714 
0.5161 

0.6115 

0.5684 

0.3317 

-0.4184 

-0.8017 

-0.7577 

0.4677 

0.6717 

0.7332 

-0.4384 

-0.6328 

-0.3598 

0.2647 


Jointed Structure 

Imaginary Part 
0.00 
-0.0082 
0.0421 
0.0728 
-0.0066 
0.0257 
0.0053 
0.0231 
-0.0613 
-0.013 
0.007 9 
0.0458 
-0.0304 

0.00 

0.3365 

-0.00796 

-0.0826 

-0.3957 

-0.4656 

0.08875 

0.0606 

0.4822 

-0.3121 

-0.4186 

-0.1236 

0.060 



Mode 3 


Mode 4 


0.2787 

0.3431 

- 0.2222 

0.1183 

0.5357 

0.5587 

0.0583 

0.3181 

0.3752 

0.2932 

0.4122 

0.0155 

-0.1589 

0.3360 

0.7448 

0.1792 

-0.4063 

-0.6329 

-0.1548 

-0.6002 

-0.8456 

0.5874 

0.5435 

-0.1980 

0.2843 


0.00 

0.1377 

-0.4765 

0.0408 

0.2585 

0.3894 

-0.1268 

0.0480 

0.1875 

0.1285 

-0.0482 

0.0035 

-0.0313 

0.00 

0.271 

0.0045 

-0.1327 

-0.3172 

-0.2363 

0.5742 

0.2825 

-0.5277 

0.2121 

0.0903 

-0.0136 


Continued 
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Mode 5 


0.3815 

0.2587 

■0.2996 

0.1225 

0.6219 

0.7381 

•0.4521 

•0.3452 

•0.1957 

•0.4472 

■0.8431 

•0.3033 

0.2891 


0.00 

-0.0224 

-0.0185 

0.0308 

0.3059 

0.4981 

0.1031 

0.1162 

-0.1541 

-0.1989 

-0.2131 

-0.1358 

0.0367 
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Table 6*10 

* Natural frequencies 
factors of Jointed 

and damping loss 
structure 

Mode 

Natural frequency 
(Hz) 

Damping loss 
factor 

1 

103.51 

0.013878 

2 

207.50 

0.000773 

3 

288.67 

0.002160 

4 

366.16 

0.009186 


5 


375.11 


0.009007 
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Fig. 6.15(b) Mode shape of jointed structre 
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substructures was done using bolts. A nut was introduced 
between the two substructures to simulate the point contact 
within the experimental limits. Upto a frequency range of 
2KHz / the beam exhibited four modes and the U— shaped struct- 
ure two modes. The FRFs taken for analysing the beam were 
for the points ^b, and \ shown in Fig. 6.16(a) and for 
the U-shaped structure the FRFs considered were for the 
points and shown in Fig. 6.16(b). 

As a first step in the coupling procedure for these 
two substructures to be coupled at points 1 and 2, the FRFs 
were taken for each substructure separately. The beam thus 
requires six FRFs and the U-shaped structure requires three 
FRFs. The modal parameters for these two substructures were 
calculated using SDF circle fit method. Table 6.11 shows 
the parameters values for the beam and table 6.12 shows 
the same for the U-shaped substructvure. Figure 6.17 and 
Figure 6.18 show the regenerated FRFs for the beam and U-shaped 
structure, respectively. The measured FRFs are also shown 
in these figures. It is clear from Figs. 6.17 and 6.18, that 
the regeneration is satisfactory and the modal parameters 
estimated are acceptable. 

The next step requires calculation of impedance as 
a function of frequency for each substructure by using (4.9). 
The receptance values for this step are taken from the 
regenerated FRF rather than the experimentally measured FRF. 
This is to avoid the noise present in the experimental FRF. 
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12.5 



(b) 



(c) 


Fig.6.16 Coupling of two substructures 




Table 6.11 : Parameters for sub-structure-l 
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Table 6.12 : Modal Parameters for the Substructure - U-shaped Structure 
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(a) FRF a* T 
b-^b 



(b) FRF a* - 
■^b^b 

Figure 6.17: Measured and regenerated FRFs for subst 'uctu-e beam 



t 



(c) FRF a* 




(.d) FRF a* , 
b^b 


Figure 6.17 
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(a) FRF a* , 
b-^b 



T? -T 17 
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Cci FRF 

( 



u 


Ficfure 6.18 
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The impedance matrix of each substructure at a particular 
frequency is obtained by inverting numerically the receptance 
matrix of each substructure. The method adopted for obtain- 
ing this inversion of a complex matrix is given in Appendix D. 


The impedance matrix for the beam is obtained by 


the 

expression given belov/ 

# 

• 





* 

•k 

* 

r * 

cs 

(o))] = 

“l 1 
-^b-^b 

* 

a 

Vb 

* 

a 



0 

a 

★ 



2^1, 

2 2, 

“ 2, 3, 



b b 

b b 

b b 



* 

* 

* 



1 

•^b^b 

01 ^ ^ 
3^2^ 

D D 

“3 3 
b b 


at frequ- 
ency 44 

(6.7) 

* 

and the impedance matrix for the U-shaped structure, [Z^(u)J 
is given by 


[ Z* (a))J =[a*(w)] 


■lu^u 


'2 1 
u n 


1 2 
u u 


'2 2 
u n 




Jat frequency 


is 


( 6 . 8 ) 

The Impedance coupling is done by using (4.11) and 
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f2,(Q)] + 




* 

KK 

b b 


* 

^1 1 

u u 

^1. 2, 

b b 

+ "*i 2 

u u 

* 


* 

* 

•k 


+ 

^2 1 
u u 

^2 2 
b^b 

" ^2 2 

u u 


'^b^b 


Z 3 2 ^3. 3, 

\^b t) b 


(6.9) 


The receptance matrix for the combined structure is 
obtained by the inversion of the impedance matrix of ( 6 , 9 ). 

[a^U)] = [Z^(a)) 3 (6.10) 


As an example, the values of receptances for the two 
substructures and the coupled structure at a frequency of 100 
Hz is listed in Table 6.13. The receptances for the coupled 
structure are obtained as point and transfer receptances for 
the location points 1, 2 and 3. Thus six FRFs of the coupled 
structxire obtained using C6.101 are the plotted in Fig. 6,19, 
along with the measured FRFs obtained experimentally. 



Ta bl e 6.13 


Values of Receptances and Impedances 
of Substructures and Coupled Structure 
at a Frequency of 100 Hz 


(i) , Substructure Beam 

Receptance ; 


FRF 

Real part 

Imaginary part 

11 

-1.134E-05 

3.081E-06 

12 

-2. 188E-06 

5.599E-06 

13 

6.066E-05 

-1.559E-06 

22 

-1.559E-05 

1.592E-06 

23 

8.137E-06 

-4.549E-06 

33 

-1.021E-04 

-4.761E-06 


Impedance : 

FRF 

Real part 

Imaginary part 

11 

-2,09E-04 

-4.393E-06 

12 

1.055E-04 

-9.429E-05 

13 

1. 198E-04 

8 .321E-06 

22 

-6.185E-04 

8.295E-06 

23 

-9.855E-06 

3.343E-05 

33 

2.556E-05 

7.572E-08 


Continued 



Table 6.13 (Continn<=.f^) » 


(ii) Substructure - U-shaped structure 
Receptance ; 


FRF 


Real part 

Imaginary part 

11 


5.322E-05 

3.812E-08 

12 


1.165E-05 

-9.542E-08 

22 


4.357E-05 

3.280E-08 

Impedance: 




FRF 


Real part 

Imaginary part 

11 


1.996E-04 

-3.64E-07 

12 


-5.335E-05 

5.747E-07 

22 


2.438E-04 

-4.541E-07 

(iii) Coupled 

Structure 


Impedance : 




FRF 


Real part 

Imaginary part 

11 


-9.434E-06 

-4.757E-06 

12 


5.214E-05 

-9.371E-05 

13 


1.198E-04 

8.321E-06 

22 


-3.747E-04 

7.841E-06 

23 


-9.855E-06 

3.343E-05 

33 


2.556E-05 

7.572E-08 


Continued 



Table 6.13 (Continued) ; 


Receptance : 


FRF 

Real part ■ 

Imaginary part 

11 

-1.263E-05 

4 . 607E-06 

12 

-1.765E-06 

1.065E-05 

13 

7.393E-05 

-1. 129E-05 

22 

-2.282E-05 

2.483E-06 

23 

6.136E-06 

-1.859E-05 

33 

1.908E-05 

1. 361E-05 
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(a) FRF |a* 

c c 


(b) |a* 2 I 

c c 



(c) 


‘ 13 ^ 
c c 


(d) 


‘2 2 
c c 


Fig. 6. 19 ; Measured and predicted FRfc's c coupled structure 
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(e) 


^2 3 
c c 


^3 3 
c c 


Figure 6,19 
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Figure 6.20 shows the theoretical FRF, 
superimposed on the experimental FRF. The theoretical results 
do not exactly coincide with the experimental results for 
the coupled structure, but the overall nature of theoretical 
FRF is acceptable. As a comparison, the theoretical and 
experimental results obtained in [ 8 ] a^e reproduced in 
Fig. 6.21. 
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CHAPTER - 7 

CONCLUSIONS 

The SDF circle fit method for modal parameter 
extraction is simple and requires little computation time. 

The present work has tried to identify the cases, where 
this SDF method does not yield acceptable results. These 
cases corresponded to system with high damping and system 
with closely placed natural frequencies. If the two closely 
placed modes are of equal strength and are highly damped 
then the circle fit method is still capable of generating 
good estimates of modal parameters. Otherwise, MDF curve fit 
method was found to be necessary to obtain the acceptable 
estimates of the modal parameters. 

For the heavily damped system, the MDF method, 
which is iterative in nature, was observed to converge 
monotonically to the correct solution. For the other case of 
a lightly damped structure with two closely placed modes 
with one mode being predominant over the other, the MDF method 
was found to be stable in only a narrow zone around the 
actual solution. This prompted the idea of reducing the 
incremental change vector ({5k} of (6.1)} by a factor e 
so that the new estimates also lie in the stable zone. 
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If sll the modal parameters are attempted to 
be modified by the MDF ciirve fit, the curve fit being done 
over the entire frequency range is very time consuming. 

Instead, it is suggested that the MDF curve fit may be 
used only for the modal parameters corresponding to closely 
placed modes and the curve fit should be done only over 
the frequency range encompassing the two modes. For all 
other modes, the circle fit method is acceptable. 

This work also attempted to verify the theory 
of impedance coupling to predict the response, of a coupled 
structure, when the responses of the individual substructures 
are known. The predicted responses, though observed to be 
similar in nature with the experimental responses, were not 
exhibiting complete informations about the modal parameters 
corresponding to some of the modes observed in experimental 
response. This may be attributed to the fact that i) the 
rotational coordinates at the coupling points were not included, 
and ii) the coupling was idealized to be a perfectly rigid 
one. With a proper model of coupling, the predicted responses 
are expected to be closer to the theoretical ones. 
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APPENDIX A 


LEAST SQUARES CIRCLE FITTING 


Let real and imaginary parts of a point on the 
frequency response function is denoted Xj, and y^ respectively. 
A number of points (say n) around a resonance is taken to 
fit a circle, using the equation given below: 


^k 


I ^k 


^k^ ^ 

■^^^k ^k^ ^k 

■2(Xk + yk^ 


(A.l) 


The summation in the above equation is for k - l...n. 
Equation CA.l) is solved to obtain a,b and c . The radius (r) 
and the centre Cx^,, y^) of the fitted circle is given by 


(a^/4) + CbV4) - c 


(A. 21 


-a/2 

-b/2 
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^pendix - B 

LEAST SQUARE ERROR MINTTMTSATTOM 


The Least square 


error e is given by 


e 


E 

over each 
frequency 
of interest 





.* 


a. 

m 


) 


^ (B.l) 


where 


_* 

cx^ are the conjugates of 


Og {measured 


cr and 
e ^ 

receptance) and (mathematical receptance) respectively. 

This error has to be minimised w.r.t. each of the 
parameters to be improved, i = l (4 n+2). 


9e 

3ir“ 

1 


- - I (a. 


* 

m 


9 a, 


m 


9k. 


+ (a 

e 


* 

m 


9 a 


m 


aki 


i = 1 (4N+2) 


0 

(B.2) 




+ iG^, then 


l^R^jk 


^1' ^3 “ ® 1 ' ^4 = "^l' ^5 = 


'^4N+2 “ ^^R^jk 
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H F. ftj, 


‘^(4N+2) ^ 


then let 


(b.3) 


(1^2 t 611^2' ^2 t 5 h , 


’''««) ^ ®(4K+2)> = 0 


(B.4) 


^ “ 1....C4N+2)) are the initial estii^ 


ates of kj 


and be obtained by SDF methods. Jk^ is small oompared to 
's for all J and is the value which are to be henceforth. 

By Taylor ‘ s theorem. 


+ ikj, + kj, 




4N+2 ap 

j=i ~9)e^ 

j 


+ higher ord^^- ^erms 


(B.5) 


By ignoring the higher order terms and usit^g Eq.(B4) 
the above equation becomes. 


4N+2 


'’I'-’^i' >^2 > + /=i 5kj k'....i, 0 


(B.61 


Using eq. (B,31, the above equation may be 


expressed 


as follows 
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3e 

3k. 


+ 


4N+2 




3k^ 3k^. 


(B.7) 


or expressed conveniently in the matrix form, 
L P J { <Sk > + {R} = i 0 ) 


(B .8) 


where any element of [p] matrix is given by 


P. . = 


3^e 


3k^3kj 


over each 
frequency 
of interest 


_* 

~ 2 a 

^ ^ d m 

^°‘e ~ “m^ 3k^3kj 


* 

3a 

m 

^ “m ^ 

^ “ e 

3 • 

3k. 

1 

m 

3a^ 

3a 

1 

3k. 3k . 

1 3 

'3kj 

‘ Tk 


_* 
a ) 
m 


(B.9) 


^tj second differentials are difficult 

tO- obtain. Also being very small, can be neglected 




3a 


m 


3 a 


m 


3a 


over each 
frequency 
of interest 


^3k 


3k. 


+ 


m 


3a 


m 


3k. 

3 


3k. 

1 


(B.IO) 
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Each element of vector {r } , 


is given by 


9e 



, * * 

2 (a - a ) 

over each ^ 

frequency 
of inter- 
est 


9 


* 


a 


m 



+ 






(B.ll) 
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APPENDIX C 


FOURIER TRANSFORMATION FROM TIME DOMAIN 


Equation (5.1) is 


x(t) = f h(T) f(t -t) dT 

^ CO 

★ 

Fourier transform X (cd) of x{t) is defined as 
X*(ui) = / x(t) dt 


Substituting (C.l) in (C.2), 


^ os -2Tra)ti " 

X (o)) = / e dt / htt) f(t -t) dT 




.By setting (t-x) as C in (C.3) 

X*C») = ” h(x) dT 7 


or 


X* (.0)1 = H* (o)) . F* (oj) 


(C.l) 

( 0 . 2 ) 


h(T)dT 

(0.3) 


dC 


(0.4) 
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APPEND ry n 

inversion of a complex matrix 

If the inverse of a complex matrix [[Aj + i[ B]]is 
given hy c j + i [ dj J , then 

([AJ+ifBj) + ([c J + i [d]) = Di;i + i[0] 

([a] [C ]- [ B ][D } + iC[B][C] + [A] [DJ)= 1 + i [ o] 

(D.l) 


From (D.l), 

CfAj [C] - [ B][D ] = fs ig (D.2) 

t[B][c] + fA][D]i = j -0 ] (D.3) 

Solving CD. 2) and (D.3I for C and D we get 

[c] = ([a] + [ B ] [ A [ B] ) “^ (D.4) 

[D] = -[Aj-^E B] C[A ]+ [B ]fA 

The real matrix inversion required above, 
was done using usual inversion process by Gauss- 
Jordan‘s method. 
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